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Dilatonic Equation of Hydrostatic Equilibrium and 
Neutron Star Structure 

S. H. Hendi * 1 ’ 2 • G. H. Bordbar 1 ’ 2 • 

B. Eslam Panah 1 • M. Najafi 1 


Abstract In this paper, we present a new hydrostatic 
equilibrium equation related to dilaton gravity. We con¬ 
sider a spherical symmetric metric to obtain the hy¬ 
drostatic equilibrium equation of stars in 4-dimensions, 
and generalize TOV equation to the case of regarding 
a dilaton field. Then, we calculate the structure prop¬ 
erties of neutron star using our obtained hydrostatic 
equilibrium equation employing the modern equations 
of state of neutron star matter derived from microscopic 
calculations. We show that the maximum mass of neu¬ 
tron star depends on the parameters of dilaton field and 
cosmological constant. In other words, by setting the 
parameters of new hydrostatic equilibrium equation, we 
calculate the maximum mass of neutron star. 

Keywords hydrostatic equilibrium of stars; dilaton 
gravity; neutron star; maximum mass 


1 Introduction 


Our observations of Supe rn ovae type la (IRiess ct al. 2004 ; 
Perlmutter et ah 1999a; Perlmutter et al. 19991 1) con- 
firm that the expansion of our Universe is currently 
undergoing a period of acceleration. But Einstein (EN) 
gravity can not explain this acceleration. In addition, 
although Einstein’s theory can explain the solar system 
phenomena successfully, when we want to study be¬ 
yond the solar system or when the gravity is so strong, 
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this theory encounters with some problems, and so 
we need to modify EN gravity. In order to improve 
EN gravity, one may add a (cosmological ) constan t 
to its Lagrangian ( Padmanabhan 2003 : Frieman 20081) . 
Moreover, we can regard other modifi cations of Ein ¬ 
stein gravity such as Lovelock gravity (Lovelock 197ll 


Lovelock 19721 Deruelle 1990; Hendi 20081 ), brane world ! 


cosmology (Demetrian 2006 : Brax 2003T Gergel^_2006 

Amarilla 2010lf. scalar-tensor theories (JordanM955 


Brans 196ll: Cai 20071 Fuiii 2003 : Sotiriou 2006 : Giddings 1993 :1 


Gregory 1997t ~i» 201 .'4 Klenac 20021: Ghodrati 20081 

and F(R) gravity (Bamba 2( if? Gognola, 20081 Corda 2009t | 
SotiBnu^.2010 • Noiiri 201ll Hendi 2012l ~ Hendi 2014al 


Momcni 2014 h . 


On the other hand, dark energy and dark matter 
have received a lot of attention in recent y ears. Theo¬ 
retical physicists introduced non baryonic (IDick 19961) 
dark matters with three models of cold, warm and hot. 
Among them cold dark matter model has the high¬ 
est agreement with the experimental observations. It 
| is worthwhile to mention that, dilaton field is one of 
the most i nteresting candidates for cold dark matter 
( Cho 199011 . In addition, one of the the best approaches 
for finding the nature of d ark energy is t aking into ac - 
count a new scalar field ( Huang 2007 : Huang 2008lf . 
Moreover, the low energy limit of string theory con¬ 
tains a dilaton field which is coupled to gravity. Phys¬ 
ical properties, thermodynamics, and thermal stabil¬ 
ity of the black object solutions in the c ontext of d ila¬ 
ton theory have been investigated before ( Tamaki 2000l 


Yamazaki 2001 : Yazadiiev 2005 : Dehghani 2008h . 


The hydrostatic equilibrium equation (HEE) plays 
crucial role in studying the evolution of the stars. This 
equation is giving an insight regarding the equilibrium 
state between internal pressure and gravitational force 
of the stars. 

It is important to note that the neutron and quark 
stars have large amount of mass concentrated in small 
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radius. Therefore, they are in the category of highly 
dense objects, the so-called compact stars. Due to 
this fact, we need to take into account the effects 
of general relativity such as the curvature of space- 
time for studying the compact stars. The first HEE 


man, OoDcnlieimer and Volkoff (TOY') (Tolman 1934; 

Tolman 19391 OoDenheimer 19391). A1 

so, the physi- 
equation have 

cal characteristics of stars using TOV 

been investigated in Refs. (Silbar 2004 

Narain 2006; 

Bordbar 2006a: Bordbar 2006bl Li 2012; 

Bordbar 20131 

Yazdizadeh 20131: Oliveira 20141). On the other hand. 


if one is interested in studying the structure and evo¬ 
lution of stars in different gravities, one should ob¬ 
tain the HEE in those gravity models. In recent 
years, the generalizations and modifications of this 
equation were of s pecial interests fo r many authors 
( Hevdarifard 2009t Astashenok 2013 : Orellana 2013 


Arbanil 2013 ]_ Doneva 201. 'll Goswami 2014 [ Lemos 2015 :| 


Yazadiiev 2014: Canozziello 2011; CaDOzzicllo 2012; 

Momeni 2015al. For more examples. 

we may note 
f(G) gravities 

that the HEE equation in f(R) and 

werediscussed in (Astashenok 2015b; 

Momeni 2015b; 

Astashenok 2015a; Abbas 20151). d-dimensional HEE 

in EN gravity was investigated in (Ponce 2000) and 

HEE of EN-A gravity with arbitrary dim 

ensions was ob- 
and higher di- 
met (GB) grav- 

tained in (Bordbar 20151). In addition, c 

mensional HEE in context of Gauss-Bor 

ity was extracted in (Zlian-Ying 2012; 

Hansrai 2015: 


B ordba r 20151) . Recently, (2 + l)-dimensional HEE was 
obtained for a static star in the presence of cosmological 
constant (IDiaz 2014T) . 

In this paper, we want to obtain modified HEE in the 
presence of dilaton field. We consider the Lagrangian 
of Einstein-dilaton gravity and a perfect fluid energy- 
momentum tensor with spherical symmetric metric to 
obtain dilatonic HEE. We also consider dilaton field as 
a correction of Einstein gravity to obtain a perturba¬ 
tive HEE. In other words, we obtain HEE for Einstein 
gravity with an additional term which the trace of small 
dilaton field. 


2 Equation of hydrostatic equilibrium with a 
dilaton field 

The action of dilaton gravity in the context of Einstein 
gravity is given by 

Ia = ± I dW^ 9 [R - 2 - V(<t>)} + I M , 

( 1 ) 

where R and $ are, respectively, the Ricci scalar and 
the dilaton field. Also V($) is the potential for 4>, and 


Im denotes the action of matter field which we consider 
a perfect fluid. Varying the action 0 with respect 
to the metric tensor g" and the dilaton field $, the 
equations of motion for this gravity can be written as 


K - 


v 2 $ 


2 d^d^- l -g^V^) 
-g^d^ + KT;, 
ldV 
4 3$’ 


( 2 ) 

(3) 


where K = In order to construct consistent so- 

lutions of the field equations 0 and (El) , we a dopt the 
approach of gravitational papers (|Chan 1995t) and as¬ 
sume that the dilaton potential contains two Liouville 
terms 


V($) = 2A 0 e 2fo<I> - 2Ae 2?$ , 


(4) 


where Ao, A, £o, and £ are constants. T his kind of 
potential was previously investigated in (lChan 19951 1. 

In the present work, we want to obtain the static so¬ 
lutions of Eq. 0. So, we assume the spacetime metric 
has the following form 


dr 2 


ds 2 = f(r)dt 2 -— — r 2 R 2 (r ) idO 2 + sin 2 ddip 2 ) , (5) 

g{r) 


where /(r), g(r) and R{r) are functions of r which 
should be determined. 

The equations 0 and © contain four unknown 
functions /(r), g(r), R{r) and 4>(r). In order to solve 
them, we consider the ansatz 


R(r) = 


— -a*(r) 


( 6 ) 


This ansatz was first introduced in ( Dehghani 2005h 
for the purpose of finding black string solutions of Ein¬ 
stein Maxwell dilaton gravity. It is notable that in the 
absence of the dilaton field (a = 0 and so R(r) = 1), 
dilaton gravity reduces to EN gravity. Also, using Eqs. 
© and © and the metric introduced in Eq. © , we 
can obtain 


<f>(r) 



(7) 


where b is an arbitrary constant and Kij = i + ja 2 . 

On the other hand, the energy-momentum tensor for 
a perfect fluid is 


T^ = (P + pc 2 ) + Pg' 11 ', 


( 8 ) 


where P and p are, respectively, pressure and density 
of the fluid which are measured by the local observer, 
and Uy, is the fluid four-velocity. Using Eq. © and 
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the metric introduced in Eq. ©, we can obtain the 
components of energy-momentum for (3+l)-dimensions 
as follows 


— nr 2 
1 o — P c 


& 


T{ = = T 3 J = -P. 


(9) 


Now, we consider the metric © and Eq. ([9]) for the 
perfect fluid to obtain the components of Eq. © with 
the following forms 


1 - r 2 AT 2 “ : 


Kpc 2 = 


+ 


crT 


KP = - 


T“ 2 r 2 ' 6 2 /C_i,i 

{g + rg'JC ip) 

(l — r 2 AT 2 “ 2 ) a 2 X 

T“ 2 r 2 6 2 /C_i,i 

g(f + rf'IC 1A ) 
r 2 //C 2 ’ 


( 10 ) 


( 11 ) 


KP = - 


(4/AT“ 2 ~2 ff /"-//'(f)') 


crT 


+ 


4/ 

(/<?)' 


a 2 5 


2r//C li i r 2 ^ 2 ^ 


where 

T = 


&\ K i,i 

r. 


( 12 ) 


(13) 


and /, g , p and P are functions of r. It is notable 
that the prime and double prime are, respectively, the 
first and second derivatives with respect to r. On the 
other hand, substituting a = 0 in Eqs. mm , one 
finds the correspon ding field equations of EN gravity 
(see ( Bordbar 2015lf for more details). 

Using Eqs. mm and after some algebraic calcu¬ 
lations, we obtain 

dP f'^p + P)^.! 2a 2 (l - r 2 AT 2 “ 2 ) 


dr 


+ 


2a 2 T 


2 / 


ATT a2 r 3 /Ci,i 


rA'6 2 /Ci,i/C_i,i rfK. 1,1 


a 2 gf" a 2 (rg' + 2 g) 

+ —— = 0.(14) 


r 3 /C4i 


In addition, one can use Eq. m to obtain /' with 
the following form 


r = 


rf(^-+KT a2 +KP) 1C 1A 

g 

, m 2 T/Ci,!/ / 


gb 2 K.—\p rICip 


(15) 


To obtain the function g(r), we consider Eq. (flOl) . 
After integration we achieve 


g{r) = 


a 2 T 


AT 0 


1 


b 2 ICiplCi t -i 
Kc 2 JC hl 


T“ 2 r 2 /Ci, : 


1 C 


Ip 2 


^2,1 


p(r, a)r Kl ' 1 dr, 


(16) 


r ^ 1,1 


where p(r,a) = ^4^4 hr which V e ff = |7r R 3 ff and 

Reff = ( ^Ic 2 3 ) r K ' 2 ’ 3 / 3ICl ’ 1 . It is interesting to note 

that in the presence of dilaton field we should replace 
r with R e ff- In other words, dilaton field can modify 
the radius of sphere into R e ff instead of usual sphere 
with radius r. It is notable that, when a = 0 (in the 
absence of the dilaton field), Reff reduces to r and also 
(as we expect) we obtain g(r) = 1— x r 2 — where 
m(r) = f 4nr 2 p(r)dr. 

Now, we consider the integral that appears in Eq. 
(ITUl) and by using of R e f /, the equation ([TUI) turns into 


9(r) = 


AT 0 


1 


b 2 JC la K ,!,_! 
Kc 2 1C i,i 


T a r 2 /Ci : 


1 C 


Ip 2 


4irr Kl 1 


M eff (r, a), 


(17) 


where we used M e ff (r, a) = J p(r,a)4'KR 2 ^dR e f f ■ It 
is notable that, M e ff ( r,a) and R e ff are, respectively, 
the effective mass and radius as results of the presence 
of the dilaton field. In obtained solution (© 71 ) and for 
consistency we use 


£o = —, 
a 


£ = «, 


An = 


6 2 /C_ i, 


Notice that A remains as a free parameter which 
plays the role of the cosmological constant. 

Now, we can obtain the HEE for dilaton gravity. For 
this purpose, we consider the Eqs. m and ED. and 
inserting them in Eq. ED- After some manipulation 
we obtain 


dP 

dr 


(Ac 2 - BJCh + C + V) 


Kr 3 flC 1/ 2 }1 lC-iplCl 1 
where A, B, C and T> are in the following forms 

A = ^ 1 //c 11 i, i^-i/ 2 ,l'lo, 2,1 M e ff (r,a ), 

3a 2 jA ' 1 , 1 , 2 ' 

2 h 2 

AT,2,4 

2 HA 


(18) 


B — — ( AK.-iplCi/ 2 ,l'y a ^2,1,2 + 


C = 2 7 “ AC^/C-r.r Ar 2 //C 1/2 ,i - 


4y 

= r 2 //C 2 ;1 /C 1/2 ,1 [Kc 2 pK?_^ +(P + pc 2 ) y] 


V 
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and also X t . h k = ^ 2 AC 2 i/" + kICijf and y + 

Krf_ir2 ]r2 
2a 2 /^1,1^-1,1' 

In next section, we continue our paper with consid¬ 
ering dilaton gravity as a correction (perturbation) of 
Einstein gravity and we will obtain the corresponding 
HEE. 


3 Dilaton gravity as a correction of Einstein 
gravity 


The interesting agreements and acceptable results of 
the EN gravity with experimental results (observations) 
guide us to consider its modification (such as dilaton 
gravity) as a correction of EN gravity. On the other 
hand, to avoid the complexity of modified gravity the¬ 
ories and obtaining credible solutions, it is logical to 
consider the dominant perturbative terms. Therefore 
in this section, we discuss the effects of small value a. 
When a is very small, we can use series expansion in 
Eq. (11811 and keep the dominant contribution term. So, 
we keep 0 (a 2 ) and ignore 0 (a 4 ) and higher orders to 
obtain 


dP 

dr 


[3 Kc 2 m(r) + 4tt r 3 (3 KP + 2A)] 
r 2 [ 3Kc 2 m(r) — 47rr (Ar 2 + 3)] ^ 

-- + o(c,‘), ( 1 . 


AC 


•i,i 


r (Ar 2 + 3) — 


3 Kc 2 m{r) 
47T 


where H is 

3 Kc 2 [l + r 2 (A + KP)] [47rA4(r) — m(r)\ 

n = --- 

87T 

9 K 2 c 6 m 3 (r) 3Kc 2 m{r)w\ 

+ 647r 3 r 4 (P + pc 2 ) 47t 

3Kc 4 m 2 (r) [12 + r 2 (2A + KP)] 

167 v 2 r 3 (P + pc 2 ) 2 ’ 


in which 


W\ — 


P [in r — 5 — A r 2 (^ — In r) ] 

2 (P + pc 2 ) 

c 2 p [(Ar 2 + l) lnr — 2Ar 2 — l] 


+ 


2 (P + pc 2 ) 


VJi = 


- (r 2 A - 1) In ^ + Kr 2 P - 1 

18Ar 2 + A 2 r 4 + 45 
3 Kr 2 (P + pc 2 ) ’ 

r 3 (c 2 p + P) [A"P (Ar 2 — 3) — 4A] In ^ 
+ 9Kc 2 r 2 p + 12 (l + Ar 2 ) 


Pr 


2A 2 


2 (Ar 2 + 3) 2 9 Kr 3 P 2 

~Kr 2 

—c 2 rp (4Ar 2 + 3) , 


A4(r) = J r 2 \n(r)p(r)dr. 


It is notable that, as one expects the first term of Eq. 
(fl9l) is the TOV equation in the presence of cosmolog¬ 
ical constant. In addition, the second term (H) is the 
leading order term of considered dilaton field as a cor¬ 
rection to EN gravity. 

We are going to continue our paper with considering 
HEE obtained in Eq. (fTTH) and obtain the properties of 
neuron stars by using of cluster expansion of the energy. 


4 Structure properties of neutron star 

4.1 Equation of state of neutron star matter 


In order to study properties of the neutron star struc¬ 
ture, we need to consider an equation of state for neu¬ 
tron star matter. The constituents of interior part of 
a neutron star are neutrons, protons, electrons and 
muons which are in charge neutrality and beta e qui- 
librium conditions (beta-stable matter) (IShapiro 1983T1 . 
In recent years, we have obtained the equation of 
state of neutron star matter using the microscopic 
constrained variational calculations bas ed on the clus¬ 
ter expansion of the energy fu nctional (IBordbar 2002 ; 
Bordbar 2004 : Bordbar 2006a ). In these studies the 


modern two- nucleon potent ials such as the new Ar- 
gonne A Vm (IWiringa 1995H and charged dependent 
Reid-93 (IStoks 199411 have been used. One of the im¬ 
portances of these calculations is the fact that it does 
not bring any free parameter into the formalism, and 
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Fig. 1 Equation of state of neutron star matter (pressure, 
P (10 35 erg/cm 3 ) versus density, e (10 15 gr/cm 3 )). 


its results show a good convergence. In this method, a 
microscopic computation of asymmetry energy is car¬ 
ried on for the asymmetric nuclear matter calculations 
which results into more accuracy in comparison with 
other semi-empirical parabolic approximation methods. 
In fact, a microscopic calculation is required with the 
modern nucleon-nucleon potentials whi ch are ex plicitly 
depend on the isospin projection (T z ) (jBordbar 1998 1. 

In this paper, for neutron star structure, we use the 
modern equation of st ate which has b een calculated 
using AV is po tential ( Bordbar 2002 : Bordbar 2004 : 
Bordbar 2006a l and investigate some physical proper¬ 
ties of neutron star structure. This equation of state of 
neutron star matter is shown in Fig. |T] 


4.2 Maximum mass of neutron star 


In order to distinguish neutron stars from black holes, it 
is important to know the maximum gravitational mass 
of a neutron star for stability against collapse into a 
black hole. In other words, it is expected that below a 
certain maximum mass, degenerate pressure due to the 
nucleons is suffi cient to preven t an object from becom¬ 
ing a black hole (IShaniro 19831) . Therefore, determining 
the maximum gravitational mass of neutron stars is of 
special importance in astrophysics. Direct ways to mea¬ 
sure the neutron star mass are observations of the X-ray 
pulsars and X-ray bursters, but there are many errors 
in these methods and so, measuring the mass of neutron 
star is not very useful and accurate. Fortunately, the 
mass of neutron stars have been deter mined with high 
accuracy us i ng the bina r y radio pulsars (IWeisberg 1984 ; 
Liang 19861 Heap 1992 : Jonker 2003 : Quaintrell 2003h . 


Here, we calculate the maximum mass of neutron 
star using the equation of state of neutron star matter 
presented in Fig. |T] by numerical integrating the HEE 
obtained in Eq. m- This leads to the neutron star 
mass and radius as a function of central mass density 
(e c ). For this purpose, by selecting a central mass den¬ 
sity (e c ), under the boundary conditions P(r = 0) = P c 
and m(r = 0) = 0, we integrate the Eq. ITUl) outwards 
to a radius r = R in which P vanishes. This yields 
the neutron star radius R and mass M = m(R). Our 
results have been given in the following figures and ta¬ 
bles. 

Our results show that the maximum mass and radius 
of neutron stars depend on the parameters of dilaton 
field (a) and cosmological constant (A). As one can 
see in Table [TJ we have obtained the maximum mass 
and radius of neutron star in the absence of cosmolog¬ 
ical constant. The results show that as a increases, 
the maximum mass of neutron star decreases. We have 
found that the maximum mass of neutron star is af¬ 
fected with variation of a for the case of a > 10 -12 . 
In other words, we can omit the effects of a for values 
lower than about 10~ 12 , where these r esults are consis - 
tent with those were obtained in Ref. ( Bordbar 2006aT) . 

Considering a = 10 -12 and A = 0, we plot the neu¬ 
tron star gravitational mass (in solar mass unit Mq) 
as a function of central mass density (e c ) in Fig. [2] 
(up diagram). Our results show that at low densities, 
the calculated neutron star mass exhibits a minimum 
(« 0.12 Mq). It can be seen that at high densities, the 
increasing of gravitational mass becomes very slow, and 
finally it approaches a limiting value (« 1.68Mq). This 
limiting value is the maximum gravitational mass of 
neutron star, and its corresponding central density is 
the highest possible value for the neutron star central 
density. A star with higher central density would be 
unstable against the gravitational collapse to a black 
hole. For mentioned quantities of a and A (a = 10 -12 
and A = 0), we also plot the gravitational mass versus 


Table 1 Maximum mass of neutron star and its corre¬ 
sponding radius for various values of a at A = 0 and 
b= 1(T 2 . 


a 

Mmax (Mq) 

R (km) 

1.33 x 10~ n 

0.22 

12.60 

1.30 x 10" 11 

0.81 

10.01 

1.20 x 10” 11 

1.37 

9.62 

1.10 x 10 -11 

1.53 

9.33 

1.00 x 10" 11 

1.60 

9.11 

5.00 x 10~ 12 

1.68 

8.55 

1.00 x 10 -12 

1.68 

8.42 
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Table 2 Maximum mass of neutron star and its corre¬ 
sponding radius for various values of a and A at b = 10~ 2 . 



a 


Mmax (ATq) 

R (km) 


A 


1.00 

X 

10" 

-12 

0.77 

6.64 

1.00 

X 

10- 

-11 

1.00 

X 

10- 

-13 

0.78 

6.65 

1.00 

X 

10- 

-11 

1.00 

X 

10" 

-14 

0.78 

6.65 

1.00 

X 

10 - 

-11 

1.00 

X 

10- 

-12 

1.55 

8.25 

1.00 

X 

10" 

-12 

1.00 

X 

10- 

-13 

1.56 

8.25 

1.00 

X 

10 - 

-12 

1.00 

X 

10- 

-11 

1.56 

9.04 

1.00 

X 

10 - 

-13 

1.00 

X 

10- 

-13 

1.67 

8.40 

1.00 

X 

10 - 

-13 

1.00 

X 

10- 

-14 

1.67 

8.40 

1.00 

X 

10 - 

-13 

1.00 

X 

10- 

-11 

1.59 

9.10 

1.00 

X 

10" 

-14 

1.00 

X 

10- 

-12 

1.68 

8.42 

1.00 

X 

10 - 

-14 

1.00 

X 

10" 

-13 

1.68 

8.42 

1.00 

X 

10 - 

-14 




R(km) 


Table 3 Maximum mass of neutron star and its corre¬ 
sponding radius for various values of a and (negative) A at 
b = 10” 2 . 


a _ M max (M 0 ) R (km) _A_ 

1.00 x lCT ia L77 8132 -1.00 x 10^ 12 

1.00 x 10" 12 L78 8133 -1.00 x 10” 12 

1.00 x 10" 11 2125 10.57 -1.00 x 10“ 12 

1.00 x 10" 13 L70 8M4 -1.00 x 10“ 13 

1.00 x KT 12 L69 8M3 -1.00 x 10” 13 

1.00 x 10" 11 L63 9T7 -1.00 x 10~ 13 


Fig. 2 Gravitational mass versus central mass density, e c 
(10 15 gr/cm 3 ), (up) and radius (down) for b = 10 -2 , A = 0 
and a = 10~ 12 . 


radius for in Fig. [2] (down diagram). Our results show 
that for the neutron star, there are a minimum gravita¬ 
tional mass (ft 0.12Mq) and a maximum gravitational 
mass ( Ri 1.68 -Mq). It is well known that in this mass 
region, the equilibrium configuration of neutron stars 
ran exist. 

On the other hand, when we consider the effects of 
ooth a and A in structure of neutron stars, the results 
ire interesting. Our calculations show that considering 
1 fixed value for A, the maximum mass of the neutron 
star increases as a decreases. Also, when a is a fixed 
value, by decreasing A, the maximum mass increases 
(see Tabic [2] for more details). Considering the effects 
of both a and A, simultaneously, we plot the gravita¬ 
tional mass versus radius in Fig. [U It is notable that 
for the special values of A and a, our results show an 
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R{km) 



Fig. 3 Gravitational mass versus central mass density, e c 
(10 15 gr/cm 3 ), (up) and radius (down) for b = 10 —2 , A = 
1CT 14 and a = 10 -11 . 


Fig. 4 Gravitational mass versus radius for b = 10 2 and 
a = 10~ 12 at A = 1.7 x 1CT 12 (up) and A = 1.7 x 1CT 13 
(down). 


abnormality in the behavior of mass versus radius for 
neutron star. In other words, for the special values of A 
and a, we encounter with an interesting behavior. For 
these cases, with increasing radius, the gravitational 
mass increases (see Fig. [4}, while for a neutron star, 
with increasing ra dius, the mass decreases, (see Ref. 
( Astashenok 2014 ) for more details). In other words, 
there are critical values for A and ct, so that, the dia¬ 
gram related to gravitational mass versus radius is not 
similar to the diagram related to ordinary neutron star. 

In order to complete our discussion and motivated 
by AdS spacetimes, we consider negative cosmological 
constant and collect the results in table [3] We conclude 
that for a fixed value of dilaton field, increasing the 
absolute value of negative cosmological constant leads 
to increasing the maximum mass of the neutron star. 


In other words, for the negative cosmological constant, 
one finds that the fraction of maximum mass per solar 
mass can be larger than 2. 

5 Closing Remarks 

In this paper, we considered a 4-dimensional spherical 
symmetric line element and extracted the hydrostatic 
equilibrium equation of stars in dilaton gravity. We 
found that for a = 0 limit, the HEE for dilaton gravity 
reduces to Einstein-A one, as one expects. Then, we 
regarded dilaton field as a correction of EN gravity and 
showed that ^ contains usual TOV equation and an 
extra dilatonic term. 
















Considering the HEE obtained in this paper and us¬ 
ing our neutron star matter equation of state, we obtain 
the structure properties of neutron star. The results 
showed that as the parameter of dilaton gravity (a) in¬ 
creases, the maximum mass of neutron star decreases. 
In other words, as a decreases, the effects of dilaton 
gravity decreases, and the results will be close to the 
results obtained in Einstein-A gravity. These results 
indicated that the dilaton field behaves as an external 
pressure and prevents the increasing mass of neutron 
star. 

In addition, we regarded negative cosmological con¬ 
stant and found that for a fixed value of dilaton 
strength, decreasing A (increasing its absolute value) 
leads to increasing the maximum mass of the neutron 
star. In other words, for both positive and negative 
values of cosmological constant, increasing A leads to 
decreasing the gravitational mass. Therefore, one may 
regard A as an external pressure which prevents the in¬ 
creasing mass of neutron star. From another point of 
view, it may be interesting to regard A as a dynamical 
pressure (jKubiznak 2012& iHcndi 2014bt) and relate it to 
the results of the present work. 

Finally, it was proved that dilaton gravity can be 
transformed to Brans-Dicke theory with a suitable con¬ 
formal transformation. In addition, it is well-known 
that a special class of Brans-Dicke theory (wbd = 0) is 
mathematically equivalent to the f(R) gravity. So, it is 
interesting to find a relation between obtained results 
of the present paper with the corresponding f(R) grav¬ 
ity and Brans-Dicke theory. We left these problems for 
future works. 
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